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Abstract 

> ■ 

We construct spherically symmetric thin-shell wormholes supported by a generalized 
Chaplygin gas in Born-Infeld electrodynamics coupled to Einstein gravity, and we analyze their 
stability under radial perturbations. For different values of the Born-Infeld parameter and the 
charge, we compare the results with those obtained in a previous work for Maxwell electrody- 
namics. The stability region in the parameter space reduces and then disappears as the value of 
the Born-Infeld parameter is modified in the sense of a larger departure from Maxwell theory. 
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<N . 1 Introduction 

> 

Traversable Lorentzian wormholes are theoretical objects which have received great attention in 
the last two decades. These objects have a throat that connects two regions of the same universe 
or two different universes [HE] and, in general relativity, they are characterized by being threaded 
by matter that violates the null energy condition [IH3] in order to allow material systems travel 
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^vq . through them. The amount of this exotic matter can be made arbitrary small [4], but at the expense 
of large pressures at the throat [5]. Traversable wormholes can be constructed [5] by using the well 
known thin-shell formalism -commonly employed to model layers in different physical contexts, 
including modern cosmology (e.g. braneworlds) and gravastars-, consisting in the cut and paste 
of two manifolds to form a new one, with a shell at the joining surface corresponding the throat, 
where the fulfillment of the flare-out condition is required. These wormholes are of particular 
interest for their simplicity, which facilitates the stability analysis, and because the presence of 
exotic matter is confined to the shell. For these reasons, they are widely studied in the literature. 
Models of wormholes with a continuous energy-stress tensor at the throat usually also require a cut 
and paste procedure to confine the exotic matter and/or obtain a suitable asymptotic behavior. 
Stability studies of spherically symmetric thin-shell wormholes, with a linearized equation of state 
at the throat, have been performed under radial perturbations( [6l-ll0| and references therein). 
Plane and cylindrical thin-shell wormholes were also considered in recent years (see, for example, 
Refs. [ME]). 
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Born-Infeld electrodynamics p3] is a non-linear theory proposed in order to avoid the infinite 
self energies of charged point particles arising in Maxwell theory, and it is the only non-linear 
theory without birefringence. The spherically symmetric solution in general relativity coupled to 
Born-Infeld electrodynamics was obtained by Hoffmann [T3]; this solution failed to be a suitable 
model for the electron, corresponding instead to a black hole. Born-Infeld type actions have 
appeared in low energy string theory |15j . leading to an increase in the interest of non-linear 
electrodynamics. Maxwell and Born-Infeld theories have electric-magnetic duality invariance [16], 
property not shared with other electromagnetic theories. The geodesic structure of Einstein-Born- 
Infeld black holes was studied in Ref. [17J. The linearized stability of spherical shells and thin-shell 
wormholes under radial perturbations was recently considered within this theory [18] . 

In the framework of general relativity, the accelerated expansion of the Universe violates the 
strong energy condition. Several models of exotic matter, proposed in cosmology |19| . have been 
also adopted in wormhole spacetimes. One of them, the Chaplygin gas [20l[21] was used as the 
exotic matter supporting wormholes [22H26] . In particular, a generalized Chaplygin gas was taken 
as the exotic matter at the throat of thin-shell wormholes in Refs. |25p 26j. 

In the present work, we construct thin-shell wormholes with a generalized Chaplygin gas in 
Einstein-Born-Infeld theory, and we analyze their stability under perturbations that conserve the 
symmetry. As in Ref. [25], we adopt the generalized Chaplygin gas equation of state at the throat 
mainly because of its interest in modern cosmology. This fluid introduces two parameters and non 
trivial complications in the equations. Here we extend to Born-Infeld electrodynamics the analysis 
performed previously using Maxwell theory [23]. We compare the new results with those obtained 
in Ref. [25]. The main outcome of our paper is that the stability region reduces in size and finally 
fades away as the Born-Infeld parameter takes values in the direction of a larger deviation from 
Maxwell electrodynamics. The paper is organized as follows: in Sec. El the wormhole construction 
is done; in Sec. El the stability of static configurations is analyzed; finally, in Sec. H] a summary is 
presented. 

2 Wormhole construction: general equations 

The action of Born-Infeld electrodynamics coupled to Einstein gravity is given by (in units such 
as G = c = 1) 

S = j d*Xy/g(-^R + L\, (1) 

where R is the scalar of curvature, g = det\g^ u \, and L depends on the electromagnetic tensor in a 
non-linear form: 

L = 4^2 ( X " ^ + \ F ^ F ° Vh2 ~ \*F m ,F*»b^ , (2) 

with F av = dfjAy — d v A a the electromagnetic tensor, *F C7I/ = \\J —ge a ^ av F a ^ the Hodge dual of 
F au , and £ a /3au the Levi-Civita symbol. The parameter b indicates how much Born-Infeld and 
Maxwell electrodynamics differ; b^ 1 is the maximum of the electric field. In the limit b — > 0, the 
Maxwell lagrangian is recovered. 

The field equations, obtained from the action (TjQ), have the vacuum spherically symmetric 
solution [THl fTT]: 

ds 2 = -iP{r)dt 2 + iP(r)- 1 dr 2 + r 2 (d9 2 + sin 2 6d(j) 2 ), (3) 
where r > is the radial coordinate, < 9 < tt and < (p < 2tt are the angular coordinates, and 
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i/j has the form: 
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with M the mass, Q the charge, and F(j, k) the elliptic integral of the first kind: F{j, k) = 
J sm7 [(l — z 2 )(l — k 2 z 2 )Y l l 2 dz = JJ(1 — k 2 sin 2 4>)~ x l 2 d<\>. The geometry is singular at r = 0; the 
position of the horizons, determined by the zeros of ip(r), have to be calculated numerically. The 
Schwarzschild metric is recovered if Q = and the Reissner-Nordstrom geometry is obtained by 
taking the limit b — > (for more details, see Ref. [17]). 

We start from the metric showed in Eq. ([3]) to construct thin-shell wormholes by using the 
Darmois-Israel formalism |27j . We need to take a radius a larger than the event horizon in order 
to avoid the presence of the singularity and the horizons. We cut two identical copies of the region 
with r > a: 

M ± = {X a = (t,r,9, V )/r>a}, (5) 



and paste them at the hypersurface 



{X/F(r) 



a = 0}, 



(6) 



to create a new geodesically complete manifold Ai = M + U M~ . If this construction satisfies the 
flare-out condition, the manifold represents a wormhole with two regions connected by a throat of 
radius a, which corresponds to the surface of minimal area. The flare-out condition is satisfied 
in our case, because ip'{a) = 2a > 0. A global radial coordinate can be defined on A4 by using 
the proper radial distance: I = ± J ^l/ip(r)dr, the signs ± correspond respectively, to ftA + 
and Mr , and the throat is located in / = 0. At the throat X we can define the coordinates 
= (T,0,(f), with r the proper time on the shell. The throat radius is a function of time: a(r). 
A Birkhoff theorem holds for the metric adopted in the construction |16| , so the geometry remains 



static outside the throat and no gravitational waves are present. Adopting the orthonormal basis 



(asin#) 1 e tp }, it is easy to obtain the second fundamental forms (or 



extrinsic curvature) associated with the two sides of the shell: 
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2v^(a) + d 2 

where a prime represents a derivative with respect to r and the dot with respect to r. With the 
following definitions: [K t .] = -fCt — K~ , K = tr[K^] = [K\] and with the surface stress-energy 



tensor S v = diag(cr, p§,Pip), where a is the surface energy density and p^, p^, are the transverse 
pressures, the Einstein equations on the shell can be reduced to Lanczos equations: 



[Ky] + Kg iS = 8ttS, 



then we have that 



and 
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It can be seen from Eq. (|10p that a < 0, which indicates the presence of exotic matter. Here we 
model this exotic matter with a generalized Chaplygin gas on the shell S. In this gas, the pressure 
has opposite sign to the energy density. Then, the equation of state at the throat can be written 
in the following form: 



where A > and < a < 1 are constants. When a = 1 the ordinary Chaplygin gas equation of 
state p = —A/ a is recovered. The dynamical evolution of the wormhole throat can be obtained by 
replacing Eqs. (JT0|) and (JTTJ) into Eq. (fT2|) : 

{[2a + 4>'(a)}a 2 + [^(a) + a 2 ]2a} [2a] a - 2^[W] Q+1 [^(a) + d 2 ] (1 ~ a)/2 = 0. (13) 

In this way, we have found a differential equation that should be satisfied by the radius of throat for 
thin-shell wormholes in Einstein-Born-Infeld theory, threaded by exotic matter with the equation 
of state of a generalized Chaplygin gas. 



3 Stability of static solutions 



In the case of static wormholes, from Eqs. (|10p and (jlip. the surface energy density and pressure 
are given by 

\AKao) 
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If static solutions exist for a given set of parameters, they should satisfy Eq. (|13p evaluated in ao: 
(2a ) a [^(a )a 2 + 2V(a )a ] - 2,4[4™ 2 ] a +VK) (1 - Q)/2 = 0. (16) 
From Eqs. (|10p and (jlip it is easy to verify the conservation equation: 

(17) 



d , n dA 

_M) +P _ = o, 



where = 47ra 2 is the area of the wormhole throat. In Eq. (|17p . the first term represents the 
internal energy change of the throat and the second the work done by the internal forces of the 
throat. We can rewrite Eq. (|17p in the form: 



which can be integrated to give 
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(19) 



This equation can be formally inverted to obtain a = a (a). Then we replace c(a) in Eq. (jlip to 
find the equation that determines completely the dynamics of the throat: 



d 2 
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(20) 
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where V(a) can be interpreted as a potential, which can be expanded in a second order Taylor 
series around the radius ao of the static solution, in order to analyze the stability: 



V(a) = V{a ) + V'(a )(a-a ) + 



V 



(°o) 



(a - a ) 2 + 0(a - a ) 3 . 



(21) 
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The first and second derivatives of V(a) are given by 



V'(a) = V'(a) + 8air 2 a(a) [a(a) + 2p(a)] 



(22) 



K"(a) = V"(a) - 8vr 2 <^ [a(a) + 2p(a)] 2 + 2a(a) [<r(a) + p(a)] 1 + 2 



ap(a) 



(23) 



where we have used that ao' = — 2(cr + p). By replacing Eqs. (fTl"|) and (fl~5|) in Eqs. (j20|) . ([22]) . and 
(f23l . we have that F(a ) = V"(a ) = 0, and 



From this last equation we obtain the stability condition for perturbations conserving the spherical 
symmetry of the geometry: the wormhole is stable if and only if V"(ao) > 0. 

By using Eq. (|16p . we can find the possible throat radii ao, for different values of the Born-Infeld 
parameter 6, the constant A, the exponent a, the mass M and the charge Q. Since the inequality 
V'iflo) > 0, which determines whether the solution with radius ao is stable, is very complicated 
from an algebraic point of view, the results (obtained numerically) are presented graphically in 
Figs. 1-6, in which standard software was used and we have chosen the most representative figures. 
The stable solutions are shown with solid lines, while the dotted lines correspond to unstable 
configurations. The regions that have no physical meaning are shaded in gray. The results present 
an important change around Q c /M, where Q c is the critical charge, corresponding to the value 
of charge from which the original metric used in the construction has no horizons. The quotient 
Qc/M only depends on the parameter b/M. For a fixed value of b/M, the event horizon has a 
radius which decreases as \Q\/M grows, and it fades out for values of \Q\/M larger than Q c /M. 
The values of Q c /M are those for which ip{rh) = and ip'(rh) = (i.e., double root of if}(r)), to be 
obtained numerically. 

From Figs. [QtoEl we see that: 

• As discussed above, for b = (then Q c /M = 1), the Born-Infeld electrodynamics reduces to 
Maxwell theory, so that the Reissner-Nordstrom solution is used in the construction of the 
wormholes, and we recover the results obtained in Ref. |25| . 

• If < b/M < 1 the behavior of the solutions is similar to what shown in Figs. [T] and [21 for 
b/M = 1 (then Q c /M = 1.02526): 

— When < a < 1 (for example, a = 0.2): 

* For < \Q\ < Q c and \Q\ not very close to Q c , there is one unstable solution for 
each value of AM a+1 , and this behavior continues when this parameter grows. 

* For \Q\ < Qc, a small range in the values of AM a+1 is found for which three solutions 
are obtained: two of them are unstable (the largest and the smallest ones) and the 
other is stable. When AM a+l grows, the behavior described above is obtained 
again: there is only one unstable solution close to the radius of the horizon of the 
original manifold. 



V"(ao) = i>"{a Q ) + 




(24) 
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Figure 1: Stability of wormholes with param- 
eters b/M = 1 and a = 0.2. In this case, 
Q c /M = 1.02526. Solid (dotted) curves rep- 
resent static stable (unstable) solutions with 
throat radius a®. Gray zones are non-physical. 




Figure 2: Same as Fig. [T] with b/M = 1 and 
a = 1. In this case, Q c /M = 1.02526. 



* For \Q\ > Q c , the range of values of AM a+1 where there are three solutions becomes 
larger; two of them are unstable (the largest and the smallest ones) and the other 
one is stable. For large values of AM a+1 there is again only one unstable solution. 

— When a = 1: 

* For < \Q\ < Q c and \Q\ not very close to Q c , we see that a bounded range of 
values of AM a+1 exists for which there is only one unstable solution. The value of 
ao/M decreases with AM a+1 until it reaches the radius of the horizon of the original 
manifold, and then no solutions are found. 

* For \Q\ < Q c , there is a range of values of AM a+1 for which two solutions are found, 
one is stable (the smallest) and the other is unstable (the largest). From a certain 
value of AM a+1 solutions are not longer present. 

* For \Q\ > Q c , there are three solutions, two of them unstable (the largest and the 
smallest ones) and the other stable. For large values of AM a+1 there is only one 
solution, which is unstable. 

Comparing the results shown in Figs. [T] and [2] with those obtained in Ref. [25], we can see 
a similar behavior in the cases 6 = and < b/M < 1, for the same values of a. The only 
difference is found when \Q\ > Q c : the smallest of the unstable solutions for < b/M < 1 is 
not present if b = 0. 

• An analogous analysis can be done for Figs. and HI corresponding to b/M = 2 (then 
Q c /M = 1.10592): 

— When < a < 1 (for example, a = 0.2): 
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Figure 3: Same as Fig. [T] with b/M = 2 and 
a = 0.2. In this case, Q c /M = 1.10592. 
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Figure 4: Same as Fig. [T] with b/M = 2 and 
a = 1. In this case, Q c /M = 1.10592. 



* For < \Q\ < Q c and \Q\ not very close to Q c , we observe that from a certain value 
of AM a+1 there is always one unstable solution. 

* For \Q\ < Q c or \Q\ > Q c , there exists a small range of values of AM a+1 where three 
solutions can be found: two of them unstable (the smallest and the largest ones) 
and the other stable. From a certain value of AM a+1 there is only one solution, 
which is unstable. 

— When a = 1: 

* For < \Q\ < Q c and \Q\ not very close to Q c , there is a bounded range of values 
of AM a+l for which only one unstable solution is present. Again, we see that clq/M 
decreases until it reaches the radius of the horizon of the original manifold, and then 
no solutions are found. 

* For \Q\ ~ Q c , there are two solutions, one stable (the smallest one) and the other 
unstable (the largest one) for a bounded range of AM a+1 , and outside from this 
range no solutions are present. 

* For \Q\ > Qc, there is always only one solution, which is unstable. 

• From Figs. [5] and El in which b/M = 5 (then Q c /M = 1.48468), we can say that: 

— When < a < 1 (for example, a = 0.2) and for any value of the charge \Q\ there is only 
one solution, which is always unstable. 

— When a = 1: 

* For < \Q\ < Q c and \Q\ not very close to Q c , a small range of values of AM a+1 
can be found for which only one unstable solution is present, with a^/M decreasing 
quickly to reach the horizon radius of the original manifold. 

* For |Q| < Q c or |Q| > Q c , there is only one solution, which is always unstable. 

This shows that for large values of b/M, the stability region disappears completely. 
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Figure 5: Same as Fig. [T] with b/M = 5 and 
a = 0.2. In this case, Q c /M = 1.48468. 
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Figure 6: Same as Fig. [T] with b/M = 5 and 
a = 1. In this case, Q c /M = 1.48468. 



Comparing these results with those obtained in the above mentioned work [25J, we observe that 
the results obtained for values of b/M > 1 differ significantly from those with b = 0. 

4 Summary 

In this paper, we have constructed spherically symmetric wormholes by using the thin-shell for- 
malism within the framework of Einstein-Born-Infeld theory, with a Chaplygin gas at the surface 
of union, where the throat is localized. We have analyzed the stability of the wormholes under 
perturbations that preserve the symmetry. The study has been done analytically and standard 
software was used to display the results graphically. We found stable solutions for a given set of 
parameters, namely the radius of the throat ao, the parameter b of Born-Infeld electrodynamics, 
the parameters A and a of the equation of state corresponding to the generalized Chaplygin gas, 
the mass M, and the charge Q. The results were compared with those obtained in a previous 
work, where a similar study was conducted using the Reissner-Nordstrom metric. For small b/M, 
there are values of the other parameters for which the solutions are stable, results similar to the 
Reissner-Nordstrom case, except that in Einstein-Born-Infeld a new unstable solution were found 
in the vicinity of the origin for large values of Q/M. As b/M increases, i.e. that the theory is dis- 
tancing itself more from Einstein-Maxwell, the stability region becomes smaller. For large values 
of b/M the stable solutions are not longer present. 
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